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Abstract 

Three topics about the application of quenched chiral perturbation theory to 
matter fields are studied. It is proved that the hairpin axial current couplings 
in quenched chiral perturbation theories do not contribute to the quenched 
chiral singularities for one chiral loop renormalization of matter field proper- 
ties. The modification of mass corrections in the chiral limit due to nonzero 
mass splittings are studied, and selection rules for hadron decays in quenched 
QCD are obtained. 
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Most lattice QCD simulations are performed under the quenched approximation, where 
the fermionic determinant is set to unity. While in the real world the u and d quarks 
have masses much less than Aqcd, lattice QCD are often simulated with m q ~ Aqcd, and 
physical results are extracted by extraplation to the chiral limit m q — > 0. As a result, 
it is important to understand the chiral limit of quenched QCD. Over the past few years, 
quenched chiral perturbation theory (Q^PT) has been formulated as the low energy effective 
theory of quenched QCD 0-^]. In Q^PT, the cancelation of internal quark loops is enforced 
by the introduction of ghost quarks (quarks with wrong - bosonic - statistics) . Since the 
contributions of ghost loops to path integrals are equal in magnitude but opposite in sign to 
those of quark loops, there will be a complete cancelation of loop contribution if there are as 
many quarks as ghosts, with the same masses. (One may also consider partially quenched 
chiral perturbation theories ||, in which the number of quarks and the number of ghosts are 
different.) Just as in ordinary chiral perturbation theory (%PT) the low energy dynamics is 
dominated by the Goldstone bosons (pions for two flavor QCD) which are quark- ant iquark 
bound states, in QxPT the low energy dynamics is dominated by the Goldstone "bosons" 
which may be quark-antiquark, ghost-antighost, quark-antighost or ghost-antiquark bound 
states. 

The main difference between xPT and QxPT lies in the the role played by the rf meson. 
In ordinary QCD, rf acquires a heavy mass due to the necklace diagrams and hence are 
integrated out in the standard formulation of %PT. Since the necklace diagrams involve 
quark loops, however, they do not contribute in the quenched theory. As a result, the rf 
have the same mass as the other Goldstone bosons (massless in the chiral limit) and cannot 
be integrated out. It does differ from the other Goldstone bosons in having an extra term 
in its propagator due the the hairpin diagram. 



where the first row/column represents the rf meson and the second row/column represents 
the rf, the ghost-antighost analog of rf, and n is both the number of quark flavors and the 
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number of ghost flavors. The first term is the usual propagator where rf and ff propagate 
independently with a single pole at p 2 = m 2 , which denotes both m v / and as it and rf 
are degenerate in quenched QCD. The second term carries a double pole and represents a 
possible mixing between rf and fj', characterized by the parameters Aq and Mq. Q 

Besides the Goldstone boson fields, one can also study the dynamics of matter fields (vec- 
tor and tensor mesons Hf3, baryons [[§HlOll, heavy mesons jTTHl^], heavy baryons [Q, etc.) 
under the framework of QxPT. In this paper, we will discuss three theoretical issue around 
the application of QxPT on matter fields. The first section will prove the (ir)relevance of 
hairpin axial couplings (which will be defined below) to quenched chiral corrections of mat- 
ter field quantities. Then the effects of mass splittings in QxPT will be discussed, followed 
by the last section which study the possibility of understanding hadron decay in Q%PT. 



I. HAIRPIN AXIAL COUPLINGS AND QUENCHED CHIRAL CORRECTIONS 

Since the matter fields are significantly heavier than Aqcd, it is customary to introduce 
the matter fields as heavy particles, i.e., static fields with definite velocities. Under the 
heavy particle formalism, the propagators of the matter fields are always (v ■ p)^ 1 , and in 
the leading order of a derivative expansion the matter fields couples to the Goldstone bosons 
through the axial currents. 

^standard = W Str [Xf, 4], (2) 

where = §(£9^ — ^d^) is the Goldstone boson axial current, "Str" denotes the su- 
pertrace over all the quark and ghost flavors, and is the axial current operator between 
the initial state Hi and the final state Hf. The form of X M depends on the matter fields in 

1 Here we are considering quenched QCD with two light flavors and we have set m u = nid- With 
three light flavors and m s ^ m q , the form of the double pole term will be more complicated. 
However, our analysis is independent of the particular form of the double pole term. 
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question. For example, the axial current of the vector meson fields V x takes the following 
form 

where g is an undetermined coupling constant, while the exact forms of the baryon axial 
current couplings can be found in Ref. |§. 

The chiral lagrangian of QxPT is complicated by the existence of hairpin axial current 
couplings, which couplings rf and ff to the flavor singlet axial current of matter fields. 

Aiairpin = ^ Str [X? * JStr^] 

^ (/)' W' 

= in 1 / a (V/)Str[^ „ M-if)- (4) 

v (/)' (») J 

After writing down the interaction lagrangians, one can study the chiral limit of different 
physical quantities by calculating the one chiral loop renormalizations. For example, the 
mass renormalization of the matter fields can be evaluated through the following integral, 

d A p 



Am 



(2vr) 4 



(coup, const.) 2 v 

(propagator of the Goldstone boson). (5) 



v ■ p 

where the sum is over all possible intermediate states. Note that the expression inside the 
square brackets are independent of p and hence can be taken out of the integral. So AM 
just depends on the form of the propagators of the Goldstone bosons. For all the flavor 
non-singlet states the propagators are (p 2 — m 2 )™ 1 , which is also the form of the "single 
pole" propagator of rj' and fj'. Putting these propagators in the integral with contribute to 
AM with the generic form 

AM ~ / — z ~ rn?. (6) 
J (27r) 4 f -pp 2 -m? y ' 

This m 3 ~ m 3 / 2 term has a singular dependence of m q which is due to low-energy chiral 
dynamics and cannot be canceled by counterterms. Since such a term also appear in normal 
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%PT, they are called standard chiral singularities. In contrast, if the contribution to the 
loop integral due to the "double pole" propagator is of the form: 

AM ~ / -^^^ViM „ Mlm. (7) 
J (2ttYv P (p 2 -m 2 ) 2 1 ' 

This term M 2 m ~ m~ x l 2 is more singular in the chiral limit than the m 3 term above and is a 
quenched artifact, which is usually called quenched chiral singularities. Since these quenched 
chiral singularities dominate over the standard ones in the chiral limit, understanding them 
may be helpful in the determination of the best extrapolation of lattice QCD results to the 
zero mass limit. 

Naively, both the standard and quenched chiral singularities may depend on the normal 
axial coupling g and the hairpin axial coupling h. It turns out that the quenched chiral 
singularities are independent of the value of h. The reasoning is very simple: the quenched 
chiral singularities come from rj' and fj' loops, and these Goldstone bosons can couple to the 
matter field through either the standard or hairpin axial current couplings. The matter 
field, made entirely of quarks and/or antiquarks, cannot couple to fj' through the standard 
axial coupling, hence 

^standard??' = ^^Qi ^standard ff = 0. (8a) 

On the other hand, rf and fj' couple at the same strength through the hairpin coupling, but 
with opposite signs. 

<?hairpin?7' = T^^h, S'hairpinTj' = —Ti'^h. (8b) 

Hence the total couplings are 

5-totaiT?' = n 1/2 (g + h), fftotai^ = -n 1/2 h. (8c) 
The one loop integral is proportional to 

^(flWal 7)' fi'total r)' +fi'total r)' fi'total fj' + fi'total fj' fi'total r)' + fi'total fj' fi'total fj 1 ) 

= n((g + hf + 2(g + h)(-h) + {-hf) = g 2 n. (9) 
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All h dependences have canceled, as claimed. 

Actually, the result can be obtained in a more elegant way. Let's define rf± = rf ± fj', 
where n' + and n'_ are orthogonal linear combinations. Notice that the hairpin axial couplings 
always couple to Str(A M ) = d^(n' — fj') = drf_. On the other hand, since the "double pole" 
propagator mixes rj' and fj', one can diagonalize it to find out its eigenstates. They turn out 
to be exactly r)' ± . 



rf rj' 



rf 



rf 



1 1 
1 1 



rj_ 



V+ V- 
( n n \ 
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(10) 



Since n'__ has eigenvalue zero, it does not propagate under the "double pole" propagator at 
all! In other words, quenched chiral singularities are solely due to rj' + loops, while the hairpin 
axial couplings couple only to drf_. As a result, the hairpin axial current couplings cannot 
contribute to the quenched chiral singularities. 

After proving that the hairpin axial current couplings do not contribute to quenched 
chiral singularities for mass renormalizations, it is time to study the possible generalizations 
and/or limitation of our result. First of all, can our analysis be extended to more than 
two flavors, with possibly different masses for different flavors? In the presence of different 
quark masses, it is more convenient to use instead of the singlet basis [rf ,ff), the flavor 
basis (mm, uu, dd, dd, ss, ss, . . .), where q is the ghost degenerate with quark q. The hairpin 
axial current couplings go through d^n^ = X^=i(?? — <?<?)? m the leading order of flavor 
symmetry breaking. On the other hand, the "double pole" propagator now takes the form 
(— A p 2 + Mq)M, where M is a In x In matrix. 
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M = 





uu 


UU 


UU 


1 p p 

1 u 1 u 


P P 


UU 


p p 

1 u 1 u 


p p 


dd 


PdPu 


PdPu 


dd 


PdPu 


PdPu 


ss 


P P 

1 s 1 u 


P P 


ss 


, PsPu 
\ 


p p 




UU 


dd 


uu 


' P P 
1 u 1 u 


P u pd 


dd 


PdPu 


p d pd 


ss 


y PsPu 


PsPd 



dd 


dd 


SS 


SS 


PuPd 


PuPd 


PuPs 


PuPs 


PuPd 


PuPd 


PuPs 


PuPs 


P, P, 


P, P, 

-r d±. d 


p, p 

rd^s 


p, p 


PdPd 


PdPd 


PdPs 


PdPs 


PsPd 


PsPd 


PsPs 


PsPs 


PsPd 


PsPd 


PsPs 


PsPs 



ss 



d^s 



qq 
qq 



qq qq 



1 1 



(11) 
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where P q = (p 2 — m 2 q ^) 1 . It is now easy to check that M. annihilate r]'_, especially if one 
notice that by the same representation r]'_ can be expressed as 

/ 1 \ 

-1 



rf. 



/1\ 
1 

vi; 



(12) 



1 

-1 
1 

v-i; 

and the last two dimensional column submatrix is annihilated by the 2x2 square submatrix 
in M.. Hence our result holds even when the quarks are not degenerate. In fact it only 
depends on the one-to-one degeneracy between the quarks and the ghosts. It is also clear in 
this representation that our result cannot be generalized to partially quenched theories, i.e., 
theories in which the number of quarks and ghosts are different. In fact it is straightforward 
to show that, in partially quenched theories with n quarks and k ghosts, the quenched chiral 
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singularities have the generic form (ng + (n — k)h) 2 , where g and h are the standard and 
hairpin axial current couplings respectively. 

So far we have been studying exclusively mass renormalization. Does our result hold 
for renormalization of other physical quantities? Before answering this question, we must 
ascertain what we mean by standard and/or quenched chiral singularities under such cir- 
cumstances. Suppose we are interested in studying the chiral correction to a certain operator 
O which acts on matter fields. The "single pole" propagator will lead to operator renormal- 
ization with the following integral form: 



d A p 



(coup, const.) 2 ^ 



P» Q^-^- 2 . (13) 



v ■ p v ■ pp 2 — m 2 



f 2 y (V(f)) {H\ fV H {{) ) 

We will refer whatever chiral correction induced by this term as standard chiral corrections. 
On the other hand, the "double pole" will induce corrections of another form: 



d A p 



(coup, const.) 2 v 

f 2 y ( H ty H (rf ( H t f y H (i)) 



Vv , Pfi -nA + nM { 



7 ? , (14) 



v ■ p v ■ p (p 2 — m 



which we will refer as quenched chiral corrections. Generically, if the standard chiral correc- 
tion takes the form of m a ~ m^ 2 in the chiral limit, the quenched chiral correction is more 
singular with the form M 2 m™" 2 ~ m^ 2 ^ 1 . It is straightforward to generalize our result to 
the following statement: Quenched one-chiral loop corrections for matter field properties are 
independent of the values of hairpin couplings. This is the main result of this section, and 
there is one important exception to this rule, which will be discussed later. 

The author does not claim to be the first to note this pattern of cancelation of hairpin 
axial coupling in quenched chiral divergences. Our result is "immediate" (to quote Steve 



2 In partially quenched chiral perturbation theories, the rj' and f( mesons do not have double 
pole in their propagators. Instead they have two different single poles, one at the original position 
p 2 = m 2 , while the other is shifted to p 2 = m? + 5m 2 . Please refer to Ref. || for details. Here 
we refer to the renormalization due to the shifted pole propagator as quenched chiral singularities, 
but in fact they are no more singular than the standard corrections unless one goes to the fully 
quenched limit. 
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Sharpe |15| ) from the quark flow diagrams in Ref. ||13|| . When one surveys the existing 
literature it is probable that the authors of previous articles on QxPT on matter field 
are aware of such a pattern (although, as far as the author is aware of, this is the first 
systematic discussion of this result). For example, Labrenz and Sharpe has calculated the 
quenched chiral singularities (Eq. (97 - 99) in Ref. [[J) for baryons in three flavor QxPT 
and indeed they are independent of the the hairpin axial couplings. So are the quenched 
chiral corrections to the octet baryon axial charges as calculated by Kim and Kim |J. Our 
rule also holds for heavy mesons, where Booth |]TTJ as well as Sharpe and Zhang [p^HU nas 
calculated the renormalizations of heavy meson decay constants, masses, bag constants, and 
Isgur-Wise form factors. In all these cases the hairpin axial couplings do not contribute to 
the quenched chiral corrections. The same can be said for Chiladze's work on heavy baryons 
[nj. For easy comparison with the literature we list the standard and hairpin axial coupling 
constants in Ref. in the following table: 



Matter Fields 

vector mesons 
tensor mesons 
baryons (N c = 3) 
baryons (N c = 3) 
baryons (N c — > oo) 
heavy mesons 
heavy mesons 
heavy baryons 



Reference 

Booth, Chiladze, Falk 
Chow, Rey 
Labrenz, Sharpe || 
Kim, Kim H 
Chow, Rey 
Booth Ml 



Sharpe and Zhang [1121 , 113 



Chiladze [14 



Standard 


Hairpin 


92 


94 


92 


9i 


D,F,H,C 




D,F,H,C 




9 


not considered 


9 


7 


9 


9' 


92, g3 


9i 



More interesting are the cases where this result actually fails. This is the case for the 
paper by Booth, Chiladze and Falk on QxPT for vector mesons ||, which is subsequently 
extended to the tensor mesons case by Chow and Rey ]/]]. For these cases the hairpin axial 
current couplings, which is denoted as in Ref. || , does contribute to the quenched chiral 
singularities. The reason for this breakdown of our result is clear. In our analysis above 
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we have assumed that the only source of quenched chiral singularities is the "double pole" 
term in the rj' and rf propagators, while for the case of vector and tensor mesons there are 
additional sources of quenched singularities, namely the possible hairpin diagrams for these 
matter fields themselves, which may lead to "double pole" propagators for these matter 
fields. However, if one explicitly setting the hairpin diagrams on these matter fields to 
be zero (i.e., if one set the parameters /io = and An = in Ref. ||), then our results 
are applicable. In other words, the hairpin axial current coupling cannot contribute to 
quenched chiral singularities which involve the rf and r( "double pole" propagators (I2 and 
I 4 in Ref. ||) but in general will contribute to the standard chiral singularities (ii) as well as 
the quenched chiral singularities which do not involve the 77' "double pole" propagator (J 3 ). 
Also note that there are more hairpin couplings for Q%PT for vector and tensor mesons, 
which are called g\ and g% in Ref. 0. The definition of g% does involve Str[A M ] ~ <9 M r/_ and 
hence our result is applicable to this g% coupling. Indeed it does not contribute to any of 
the quenched chiral singularities. On the other hand, the definition of g\ does not involve 
Str[A M ] and hence is not constrained by our result. 

In summary, we have proved that the quenched chiral corrections to matter field prop- 
erties are independent of the values of the hairpin axial current couplings. Besides serving 
as a useful check of the correctness of Qx?T calculation results, this observation may be 
useful in the eventual determination of parameters in the QxPT lagrangians. All the cou- 
pling constants appearing in the QxPT lagrangians are nonperturbative quantities which 
cannot be calculated from first principles from quenched QCD and has to be extracted from 
lattice results (just as the paramaters in standard xPT lagrangians have to be extracted 
from experimental data). While one can extract the "double pole" parameters from lattice 
results on Goldstone boson systems, it may be difficult to disentangle the effects of the 
standard axial couplings from their hairpin counterparts. Our result suggests the procedure 
that one should extract the values of the standard axial couplings from the quenched chiral 
corrections, and with those pieces of information one can proceed to extract the values of 
the hairpin axial couplings from the standard chiral divergences. 
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II. MASS SPLITTING AND CHIRAL MASS CORRECTIONS 



Chiral perturbation theory, both the standard and quenched versions, provides a system- 
atic framework to describe low energy interactions between hadrons. Before one constructs 
a chiral lagrangian, however, one must decide which fields to include. There is no clear 
cut criteria on which fields should be included, but the conventional wisdom suggested that 
one should at least include all fields which are related by light flavor SU(n) (n = 2 or 3, 
depending on the scenario), as well as all the fields which are degenerate in the large N c 
limit. If the hadron in question contains a heavy quark, then we should also include all 
the fields related by heavy quark spin symmetry as well. For example, in Ref. || both the 
vector meson octet and the singlet appeared in the chiral lagrangian, as the octet and the 
singlet are degenerate in the large N c limit. In fact they can be conveniently combined into 
a U(3) nonet. Similarly in Ref. |J the baryon chiral lagrangian contains both the octet and 
the decuplet, which are again degenerate in the large N c limit. It is well known that, in 
standard ;\PT, the failure to include the baryon decuplet will lead to huge chiral corrections 
to baryon octet, while the chiral loops in a with both the octet and the decuplet is 
much smaller [|16[]. Physically these degenerate states are strongly coupled to each other, 
and the failure to include them will lead to an incomplete representation of the interaction, 
causing a badly-behaved loop expansion. 

In both the real world and its quenched approximation on the lattice, however, all these 
symmetries (chiral, heavy quark, and large N c ) are broken and the degeneracies are only 
approximate. In (Q)%PT, one often (though not always) treats these mass splittings as small 
perturbations and include them leading to the leading order. While such approximations are 
usually justified, they are not entirely satisfactory as it is not clear if the higher order effects 
will change the chiral behavior drastically. The question can be put more impendingly in 
the following way. Let's take the chiral corrections to hadron masses as an example. Since 
the one loop diagram has two axial current couplings, each giving a factor of 1//, one can 
obtain from dimensional analysis that 
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AM ~ nf/f\ (15) 

simply on the basis that the Goldstone boson mass m is the only other mass scale in the 
theory. [] In other words, the chiral behavior of the mass correction cannot be changed 
unless some other mass scale is introduced in the theory. In QxPT, the new mass scale is 
Mq , and the quenched chiral singularities are proportional to M^m, which dominate in the 
chiral limit. Now the mass splitting A between states is another new mass scale and can 
potentially change the chiral behavior. The question is, will it? 

Fortunately, it is possible to study the Feynman integral exactly. This has been done in 
%PT in Ref. [|l6j for baryons (though the treatment is general and can be applied to other 
systems as well), and the generalization to QxPT is trivial. 

Let's review the analysis in Ref. [16|, which study the integral f] 



dp p 1 



(16) 



f 2 J (27r) 4 v ■ p — A p 2 — m 2 ' 

where we will assume A > so that the intermediate state is heavier than the external state 
and hence the integral does not have any cuts corresponding to real particle decays (which 
will be the subject matter of the following section). Then the integral can be evaluated 
exactly. It is instructive to separate the result into two parts: 

I = + (17) 

The analytic contribution 

/ (a) = -7^(f A 3 + ^ Am 2 ) (18) 
±2ttj z 671 71 



3 Note that chiral logarithms like m z In appear only when Z is an integer, so that the effect of 
shifting the subtraction point u can be canceled by a counterterm analytic in m q ~ m 2 . In this 
case Z = 3, and chiral logarithms are forbidden. 

4 There is a factor of i difference between our definition of I and that in Ref. 



12 



can be canceled by counterterms, while the nonanalytic contribution cannot be modified by 
counterterms and is genuine signature of low energy chiral dynamics. 

i f^(m 2 - A 2 ) 3 /2 gee- 1 f + A( A 2 - |m 2 )i In 4, A < m; 

jM = v ^ ~ ( 19 ) 

12tt/ 2 | _1 (A 2 _ m2) 3/2 ]n A-V|gg + A(A 2 _ | m2) l ln ^ A > 771, 

where /z is the subtraction point. This nonanalytic contribution is plotted in Fig. 1, with [i 
set to be equal to A and 2 A. It is clear that, while /( na ) does behave like m 3 at large m for 
both values of /i, 



j(na) 



2 



o 3 2 . . m 

m m A In — - + . . . 

2tt /i 2 



m > A, (20) 



127T/ 2 

there are nontrivial corrections in the chiral limit due to a non- vanishing A for [i = A, while 
the chiral behavior is much better-behaved with /i = 2 A. It should be understood that 
changes in [i generate extra analytic terms (m 2 and m° terms in this case), which should 
be absorbed into To see it more explicitly one can perform a Taylor expansion around 



m = 16 



j(na) = 1 1 



-A 3 + -Am 2 n -^r -In — — + O — 

2 ; 4A 2 8 A 4A 2 V A 2 



12vr/ 2 tt 

The ln/^ 2 /4A 2 term is actually analytic in m 2 which can be absorbed into the analytic 
part and canceled by counterterms. It vanishes, however, at fi = 2A, explaining why the 
/i = 2A curve approaches the chiral limit so nicely. The leading nonanalytic contribution 
is m 4 lnm 2 ~ m 2 lnm q , quite different from the m 3 behavior when A = 0. In fact the 1/A 
suppression factor is reminding us that this is a higher dimensional operator produced on 
integrating out the heavier intermediate state. 
The corresponding integral in Q%PT is 

%_ r d 4 p p 2 -A oP 2 + M 2 
q f 2 J (2tt) 4 vp-A (p 2 -m 2 ) 2 ' 1 ' 

where the subscript q stands for "quenched". As discussed above, the quenched chiral 
singularities comes from the Mq term. So for simplicity we will set Aq = is the discussion 
below, and 
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_ iM 2 r d 4 p p 2 1 _ 2 dl 

q ~~P I (2nYvp- A{p 2 -m 2 ) 2 ~ °d^ { ' 



We are only interested in the nonanalytic part I^ na \ which is 
M 2 

P I _3A l n m!^ ACA 2 _ ^2^1/2 1^ . 

A+VA 2 -m 2 ' 

We see the expected / ~ m ~ m^ 2 nonlinear behavior for large m. 



12 J r/0_|| ln ^ + J. (A2 _ m 2 ) ,/2 ln |^g i A > m _ 



j(na) 



1 



127T/ 2 

And for small m, 



2 



3 3 T71 

-M 2 m- — M 2 A In— + 

2 Z7T /i z 



m > A. (25) 



J(na) _ < 



127T/ 2 



3A , /i 2 3 m 2 , m 2 ^ . 1 , 
2^ ln 4A^-47X ln 4A^ + 0( A^ ) 



(26) 



so the leading nonanalytic term is m 2 \nm 2 ~ m q \nm q , again quite different from the 
m ~ mj/ 2 behavior with vanishing A. Another interesting point is shown in Fig. 2, where 
1^ is plotted with /i — A and 2 A. In this case a change in /x just shifts the curve vertically 
but do not change the "shape" at all, as dictated by the Taylor expansion above, where the 
coefficient in front of the ln/i 2 term is just a constant in m. Note that such a vertical shift 
would not affect extrapolations to the chiral limit at all as the "shapes" of the curves are 
identical. 

What does the analysis above imply for real lattice QCD simulations? Note that the 
nonlinear behavior in Fig. 2 appears around m/A ~ 0.5. So the implication of Fig. 2 is 
twofold. Firstly, it suggests that, for m > A/2, the extrapolation of hadron masses to the 
chiral limit should be linear, a fact which is consistent with existing lattice results. On the 
other hand, it predicts the breakdown of such linear behaviors at m < A/2, a regime which 
is not yet well-explored. In real lattice QCD simulation, the application of this analysis 
is complicated by the existence of more than one possible intermediate states, and hence 
more than one A. However, one can still make some qualitative predictions. For example, 
since the N-A splitting is around 300 MeV, one expects the extrapolation for the nucleon 
mass departs from a linear behavior at around 150 MeV (about the physical pion mass). 
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On the other hand, the small p-u splitting (17 MeV) means that probably a simple linear 
extrapolation is sufficient in the vector meson sector. Certainly more detailed studies along 
this line, especially for physical quantities other than masses, will be very useful to the 
extraction of physical information from lattice data. 

III. HADRON DECAY IN QUENCHED QCD 

In this previous section, we have studied two point functions with the intermediate state 
Xi heavier than the external state X e , i.e., A = Xj — X e > 0. In this section, we will 
consider the opposite, where A < and the intermediate state is lighter than the external 
state. In particular, we are interested in the case which A < —m, and at large Euclidean 
time the two point function will be dominated the Xj7r state, resulting in a cut in the two 
point function. However, often these decays involve creating a quark-ant iquark pair from 
the vacuum. The question is, can hadron decay be studied on lattice under the quenched 
approximation? 

In this section, we are going to derive "selection rules" for the hadron decays under 
quenched approximation. For example, an obvious selection rule is the daughter hadrons 
cannot contain any valence quark which does not appear in the mother hadron as valence 
quark, as these new valence flavors have to from an internal quark loop. We will see below 
that quenching does impose strong constraints on the possible decay mechanisms. We will 
focus at system with the standard and hairpin axial current couplings (denoted by g and h 
respectively), but no other kind of couplings (i.e., g\ and gs in Ref. are out). Then there 
are two classes of diagrams representing hadron decays: 

(1) Diagrams which all the vertices are standard couplings g, and no hairpin propagators. 
The corresponding integrals will be proportional to g 2 . Since these decays have clear analogs 
in standard (unquenched) QCD, we will call them "standard decays". 

(2) Decays with one hairpin axial vertex (integral ~ gh), and/or diagrams with the 
"double pole" propagator on the Goldstone boson line (integral ~ g 2 M^). These decays may 
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either have analogs in QCD or be merely quenched artifacts. We will call these "hairpin 
decays" . 

We will first consider standard decays for meson. It is trivial to see that OZI suppressed 
decays, like <fi — > pir, cannot take place under the quenched approximation. However, not 
all OZI allowed decays can survive quenching. Decays through Fig. 3a, which contains an 
internal quark loop, are not allowed in quenched QCD, while decays through Fig. 3b are 
preserved. Note, however, that the quark and the antiquark from an external meson line 
in Fig. 3b annihilate, demanding that the external meson must be flavor neutral. However, 
isospin symmetry, which is a good symmetry of Q%PT, can rotate flavor neutral states into 
flavor states unless the flavor neutral state has isospin zero. Hence, we have the selection rule 
that for the only possible standard meson decays are OZI allowed decays of I = mesons. 
The amplitude of such decays in quenched QCD is identical to its counterpart in standard 
QCD. An example of such decay is uj — > pn in the vector meson sector, and its tensor meson 
analog f 2 — > a 2 n; unfortunately these processes do not happen in the real world where the 
total mass of the daughter particles is larger than the mother particle. It turns out that 
light mesons, which form U(3) nonets, usually mix ideally so that the 1 = mass eigenstates 
are the ss state and the uu — dd states. The ss states, like the (f> vector meson and the f 2 
tensor meson, can either decay through OZI suppressed modes like — > p7r, which is clearly 
disallowed in quenched QCD, or modes like — > 2K, which is described by Fig. 3a and 
hence also do not survive quenching. And while the uu — dd state, like the ui vector meson 
and the f 2 tensor meson in principle can decay through Fig. 3b, they are so light that they 
cannot decay into the 1=1 state in the same nonet (p and a 2 , respectively). So the only 
possibility for standard decays in the meson sector are decays of an / = state from a highly 
excited nonet to two lighter 1=1 states in other nonets. These decays cannot be described 
by QxPT as they involve transitions between different nonets, and it is probably difficult to 
really observe these decays on the lattice as these highly excited states may be very dirty. 
In passing we note that heavy mesons do not have standard decays in quenched QCD. 

In the baryon sector we are interested in studying decays like A — > Ntt and — > A c tt. 
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The former decay has been studied in detail in Ref. ||, from which we will borrow heavily. 
The main complication in studying these baryon decays is that the initial and final baryons 
have different isospin wave functions to satisfy the spin-statistics theorem. For concreteness, 
Ref. H has chosen to study the decay of A ++ , which in standard QCD the only possible 
decay mode is: 

A ++ = uuu — > ud + duu = 7r + + p + (27) 

Note that a dd pair is created. So naively A ++ decay is prohibited by quenching. This 
argument, however, was shown to be wrong in Ref. 

Standard decays for baryons can go through either Fig. 4a, which has an internal quark 
loop and hence does not survive quenching, and Fig. 4b, which does. We will denote the 
amplitude of A ++ decaying through Fig. 4a as T q , where q is the quark on the internal loop, 
while T is the amplitude of its decay through Fig. 4b. In the unquenched world, the total 
decay amplitude is 

punquenchcd = Tq + g = T + T u + T d . (28) 

Q 

On the other hand, fermionic statistics of the quark forbids the existence of an intermediate 
uuu octet. Hence we have 

r = -r u , (29) 

and hence 

punquenched = ^ = p ? (3^ 

The A ++ decay amplitude after quenching is 

r quenchcd = ^ + £ ^ = p Q + ^ + r _ + + p. (3^ 

The cancelation of quark loops and ghost loops demands that 

r, + r 9 - = 0, (32) 
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which gives 

pqucnchcd p^ p punquenched (33) 

So the decay amplitude changes sign under quenching. A more illuminating way to present 
the result is to note that 

pquenched = ( p Q + p^ + p_ + ^ + p _) = p_ = p^ (34) 

as the sum in each pair of parentheses vanishes. So in quenched QCD A ++ decays to 
(uu) (uuu) , which clearly has the opposite sign as the QCD decay. 

How does this lesson generalize to other baryon decays? We will focus our attention 
to the orbitally unexcited baryons, which are the states with most physical interest. In 
other words, we will consider the octet-decuplet light baryon system, and the antitriplet- 
sextet heavy baryon system. Now we make the phenomeno logical observation that, while 
baryon decays through pion emission are possible in both of these systems (e.g., A — > Nn 
and E£*) — > A c 7r), in neither of these system can baryon decay through K or i] decay as 
physically these mesons are too massive. We will make the assumption that this pattern 
continue to hold in quenched QCD. Then Fig. 4b demands that, for decay through pion 
emission to take place, the initial state must has at least two light valence quarks. So 
the selection rule is: standard decays for baryons are possible if the decaying baryon has 
at least two light valence quarks, and the list includes A — * N and — > A( c ). And for 
all of these decays, the quenched amplitude has the same magnitude but opposite sign as 
the unquenched counterpart. On the other hand, orbitally unexcited S(q)'s and ^(q)'s do 
not have standard decays under the quenched approximation. Our rule does not apply to 
orbitally excited baryons, which may have enough phase space to decay through K and rj 
emission, but again studying these states on the lattice may be quite difficult. 

After studying the standard decays, let's now move on to the hairpin decays. Hairpin 
decays either have a "double pole" Goldstone propagator, leading to an integral proportional 
to Mq g 2 (the hairpin axial coupling h cannot appear; see Section 1), or are interference effects 



between the standard coupling g and the hairpin coupling h, i.e., proportional to gh. Note 
that the h 2 term is absent here as in that case the Goldstone boson line would form a quark 
loop. In other words, the hairpin vertex h can cross interfere with g, but not self interfere, 
which is an evidence for non-unitarity. 

The important observation here is, hairpin decays always have the form X e — > Xtf', as 
only i]' couples through hairpin coupling, or propagate through the "double pole" propagator. 
Since rj' is a singlet under flavor SU(2) and SU(3), however, it means that hairpin decays 
can take place only between multiplets with the same representation under the flavor group. 

In the view of this result, the most natural places to look for these hairpin decays are 
transitions between a heavy quark doublet, like D* — > Drf and S* — > X c 7/. Since the two 
states in a heavy quark doublet is related by just flipping the heavy quark spin, they always 
have identical flavor representation. As discussed above, there are no standard decay for 
D* , and this hairpin decay is the only possible mode under the quenched approximation. 
The £* baryon can decay to A c 7r, which would probably dominate over this hairpin mode, 
as the latter in suppressed in the large N c limit |TIJ . In the real world, moreover, £* is only 
80 MeV above S c , so probably this hairpin decay cannot take place after all. 

The strange quark is not massive enough to be considered a heavy quark, and while the 
"strange quark doublet" S* and E have the same representation under flavor SU(2) (both 
have 1 = 1), their flavor SU(3) representations are different (S is in the octet, £* is in 
the decuplet). Hence the hairpin decay S* — > £r/ is flavor SU(3) suppressed. Lastly, note 
that hairpin decays from the tensor meson nonet to the vector meson nonet are allowed 
under this selection rule, and in principle one can study these decays on the lattice. These 
decays, however, are not described within the framework of (Q)xPT, where the tensor meson 
number and vector meson number are separately conserved. 

In conclusion, we have chosen that quenching provides severe "selection rules" on hadron 
decays. Studying decay on lattice are important as it provide one of the most direct way 
to determine the values of the coupling constants in the QxPT lagrangian (just as exper- 
imental decay rate can lead to determinations of coupling constants in the standard chiral 
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lagrangian). One can obtained C ~ ^antt from A — >N7r, and g% ~ 5 , s C*)a 0] from 
£(*) y\ c7T _ Regretfully, hadron decays are not observable in most of the lattice QCD 
simulations in the present generation as the light quark masses are still relatively large, and 
hence m > — m-N- But it is probable that, with the improvement in algorithms and 
increase in available computational power, probing the small m region would be a feasible 
endeavor in the not-too-distant future. 
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FIGURES 

FIG. 1. Plot of -12vr/ 2 /( na )/A 3 vs m/A. The dashed curve is obtained with 11 = A, while the 
solid curve is obtained with [x = 2 A. 

FIG. 2. Plot of -127r/ 2 / (na ) / (M(jA) vs m/A. The dashed curve is obtained with fx = A, while 
the solid curve is obtained with li = 2 A. 

FIG. 3. Diagrams contributing to two point functions of a decaying meson. Fig. 3a contains an 
internal quark loop and does not survive quenching, while Fig. 3b does not contain internal quark 
loop and contributes in quenched QCD. 

FIG. 4. Diagrams contributing to two point functions of a decaying baryon. Fig. 4a contains 
an internal quark loop and does not survive quenching, while Fig. 4b does not contain internal 
quark loop and contributes in quenched QCD. 
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